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Abstract. We define modulated replica symmetry breaking (RSB) schemes which combine tree- and wave- 
like structures. A modulated scheme and unmodulated RSB are evaluated at 1-step level for a semicon- 
ductor model with antiferromagnetic Korenblit-Shender interaction. By comparison of the free energies we 
find evidence that a T = phase transition in the ferrimagnetic phase leads to a transition between the 
different RSB-schemes. An embedding factor of Parisi block matrices with sublattice-asymmetrical size is 
employed as a new variational parameter in the modulated scheme. 

PACS. 68.35.Rh Phase transitions and critical phenomena - 75.I0.Nr Spin glass and other random models 



1 introduction 

Hierarchical tree-structures and replica symmetry break- 
ing (RSB) are celebrated features in the theory of 
magnetic systems with random and frustrated interactions 
of infinite range. Their role for short-range spin glasses 
was hotly debated recently . A puzzling question over 
many years concerned the existence of variants or alterna- 
tives for the Parisi i?5i3-scheme. The latter proved to be 
very robust however; it is hard to find relevant variables 
which perturb or change the scheme. 

In this article we define and apply i^S* i?-schemes which 
show a wave-like modulation in addition to the tree-struc- 
ture We work in the context of a two sublattice 
infinite-range interaction model, defined and analyzed in 
replica-symmetric (i?5')-approximation by Korenblit and 
Shender [51 (ifS'-model). The ifS'-model successfully of- 
fered the description of transitions between spin glass and 
antiferromagnetic order (or ferrimagnetic order in a field) 
in spite of an infinite-range interaction. Spatially stag- 
gered order is allowed by letting the interaction work only 
between different sublattices. The field of application in- 
cludes two-component magnets as well as standard anti- 
ferromagnets, where staggered magnetic order defines sub- 
lattices. Range-free interactions cannot distinguish spatial 
positions and consequently unite the mean field picture of 
sublattice systems with another class of systems having 
just an equal number of mutually interacting but arbi- 
trarily placed A- and B-spins. The KS-modc\ effectively 
mirrors the phase diagram of the SK-mode\ ^ to the anti- 
ferromagnetic side, still allowing to retrieve ferromagnetic 
solutions. Intra-lattice interactions are a less relevant de- 
tail 0; they can yet be included and dealt with in a refined 



KS-mode\. Transitions from spin glass to ferrimagnetic or- 
der, driven by antiferromagnetic interactions, are frequent 
physical phenomena and concern a wide range of different 
microscopic models. Even in cases when quenched disorder 
is weak or absent, spin glass models can have the power to 
mimic behaviour of clean but geometrically frustrated sys- 
tems jTliniin] ■ Beyond the present application to antiferro- 
magnetic instabilities and ferrimagnetic glassy phases, our 
results suggest that the modulated i?S'i?-schemes could 
also serve as a basis to describe excited states in finite- 
range spin glasses. 



2 The two-sublattice spin glass model with 
competing antiferromagnetic- and 
ferrimagnetic order 

A class of Hamiltonians, for which the modulated RSB- 
scheme is constructed, is given by the KS-uiodel in an 
external field h{r) 
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where the partially frustrated random interaction JiA.jB 
is chosen to obey a Gaussian distribution P{JiA,jB) = 
exp{-N{JiA,jB + Jaf/Nf/{2J^))^N/{2TiJ). Centered at 
a negative mean coupling {Jij^j^) = —Jaf < 0, the sub- 
lattice-interaction permits glassy antiferromagnetic order. 
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Glassy ferrimagnetic order with lifted A B symmetry 
results for example when a homogeneous field is applied 
or when spins of different lengths (different spin quan- 
tum numbers in quantum models) happen to be located 
on different sublattices. Modulated RSB should also be 
considered for model classes including interactions of dif- 
ferent types of localized spins, for example tight-binding 
electron spins coupled to ionic spins 5; initially mobile 
carriers, which localize due to their interaction with ionic 
spins, may not be able to fit the RSB glassy order. This 
example reaches far beyond the classical model (1). In 
the present work we focus exclusively on model (1) with 
minimal inequivalence of sites (such as being of A- and 
-B-type), which requires a hybrid modulated form of RSB 
in a solvable classical model and hence reveals a coupling 
of replica- and real space. 

In replica theory |10| . which we use here, all spin vari- 
ables acquire a replica- index a, S — *■ S"' . After elimina- 
tion of the microscopic spins the corresponding effective 
Lagrangian of the replica theory jHj is given in terms of 
Hubbard-Stratonovich fields jS]. The SK-m.ode\ interac- 
tion requires one such field [HI, Q""'^. Its statistical average 
Qab = ]^Qab^^ ^ {Sf S'^) describes glassy order [HE] in ad- 
dition to a homogeneous magnetization M = (S'f), which 
can be finite in case of partial frustration. The KS-mode\ 
however involves for each sublattice k = A, B a, magneti- 
zation = {SfJ and = {Q"/), and a field Q'^''' 
which couples the sublattices The averaged matrix 
Qab = —i Qs = —i{Q^) turns out to be equal to Qa+Qb, 
where Qa and Qb inevitably show sublattice-splitting of 
their entries qA 7^ qB in ferrimagnetic phases, together 
with \Ma\ 7^ lAfsj. The size of their block-diagonal ma- 
trices, characterized in RSB by a Parisi parameter m 
may also develop a sublattice-asymmetry. Thus, at 
\RSB-\eve\ two order parameters for each sublattice A or 
-B, hence {qiA, q2A, qiB, 92b}, denoting matrix elements of 
Qa and Qb, and two Parisi-parameters m^,TOB need to 
be considered. 

A simple illustration for the matrix Qab = —i Q3 is 
displayed in Fig^ Note that for simplicity it is not shown 
that elements on the diagonal vanish (while, for example 
in fermionic spin glasses, these elements equal 1 at half- 
filling; this detail can trivially be accounted for in the 
trace formulas below, but is of no relevance for our present 
application). 

Unmodulated replica- symmetry breaking scheme: In the 
unmodulated IRSB scheme one chooses tua = ms- 



3 Modulated Replica Symmetry Breaking 
Schemes 

We consider the superposition of two diagonal Parisi-block 
matrices with sizes tua x tua and tub x ruB such that a 
rational embedding factor 7 denotes the number of smaller 
blocks (let ruA < tub) fully embedded inside larger m^- 
boxes, i.e. without being intersected by the latter ones, 
divided (normalized) by the number of larger boxes inside 
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Fig. 1. Simple examples of modulated saddle point matrices 
Qab = Qa + Qb with mA = 3, = 5, and n = 15 (left 
Figure) and = 7, = 11, and n = 77 (right Figure) in 
one-step RSB. Four different regions are shown with entries 
qiA + qiB, generated by the overlap-areas of niA x niA- and 
mB X m_B-sized block matrices along the diagonal, with qiA + 
q2B and q2A + 91s generated by their nonoverlap-regions, and 
q2A -\- q2B belonging to the superposition of off-diagonal A— 
and B-elements. 

the entire (n x n) host matrix (the example of Fig^show 
7=1 (left figure) and 7 = f (right figure)). 

We distinguish now a single embedding scheme, for 
which the host size n is restricted to be the least common 
multiple of niA and ms, and multiple embedding schemes, 
which align /c-times such structures along the diagonal. 

For all of n,mA,mB,n/mA, and n/mB integral and 
mA < mB, the embedding parameter 7 is given by 

^ {n/mA~ (n/mB - l))/{n/mB)- (2) 
Let us now turn to the important replica limit n — > 0. 



3.1 Single-embedding modulated scheme (SMS) 

When the limit n — > is approached, as required in replica 
theory, Eq. (O and the integral constraints associated with 
it must be relaxed. This allows to obtain a finite nontriv- 
ial free energy. To this end we choose 7 as a free vari- 
ational parameter like mA and m^. To obtain the free 
energy F for this modulated i?S'i?-scheme we need the 
trace of the square of such super-imposed Parisi matrices 
(see EqEl below). The free energy involves the limit as n 
goes to zero of (1/n) times a sum of traces over the Q^- 
values. As for standard Parisi matrices [T|P]. tr{Ql} = 
n(TO„ — l)qi^ + n{n — TOk)'Z2ki k = A,B, so for a physi- 
cal finite free energy one requires that trjQa}^ also scales 
like n for small n. A set of self-consistent equations is de- 
rived by the condition that m^, ms and 7 extremize the 
free energy in the replica limit. The idea is thus to find 
a function of mA, rriB and 7, which agrees with the trace 
of Q\b f*-"" integer- valued parameters mA,mB,n allowed 
by the scheme, satisfies the integral constraints ahead of 
Eq.©, and whose analytical continuation for small n is 
linear in n. This goal is achieved by evaluating all overlap 
contributions to tr{Q^} and reexpressing the results for 
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each of the four different overlap regions in terms of the 
embedding factor 7. 

The trace formula for for arbitrary block-sizes uia < 
niB, which all nest the host-matrix of size n, for the single- 
embedding scheme, is expressible in terms of 7 as 

(7) = tr{Q\g} = n{mA - l){qiA + Qisf 
+n{mB - mA){q2A + qisf + "-("- - ■mB){q2A + 925)^ 
+2{qiA- q2A){qiB - q2B)(l)h), (3) 

where the dependence on the embedding factor 7 is con- 
tained in 

(t>il) = ^ ■ri)imB - (1 + l)mA) 

[{niB - 2n)mB + mA{2n{l 7) - (4 + -f)mB))]}. (4) 

The result represented by Eq. (|3Jl , together with (f> given 
by Eq.Q, holds for all integral and non- integral (rational) 
values of 7 allowed by the construction. Let us consider 
a few examples, using this division into two classes with 
either integral or non-integral embedding factors 7: 

i) there exists a subset of matrices, where each tub- 
block hosts the same (integral) number of m^-blocks. The 
left hand side of FigQ] shows one example with {niA — 
3,mB = 5,n = 15,7 = !)■ Further examples of this class 
are (2,7,14,3), (4,7,28,1), (3,8,24,2); 

ii) the right hand side of FigQ] presents one example 
for the matrix-class having non-integral embedding factors 
with {ttia = 7,mB — 11, n — 77,7 — 5/7). Let us add 
further examples by (3,4, 12,2/3) and (7, 17, 119, 11/7). 

One can see that S'-contributions from all overlap re- 
gions can be expressed in terms of m^, m^, 7, while 
the number of these overlap regions depends explicitly on 
the host matrix size n. This feature guarantees the finite 
replica limit of the free energy. 

3.2 Multiple-embedding modulated scheme (MMS) 
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Fig. 2. Example of a modulated saddle point matrix belonging 
to the multiple-embedding scheme: {Q3) /i = {Qa) + {Qb) with 
niA = 4, TUB = 6, and n = 36 in one-step RSB. 

7 is thus restricted to values 7 > — 1 in contrast to the 
SMS. The MMS also differs essentially from the SMS 
by the fact that the (A,B)-symmetric limit mA = = 
m reduces it to the unmodulated li?S'i?-scheme. Another 
possible variant of MMS, where the number k = n/ni of 
repeated S'M S'-structures is varied, is discarded, since k 
cannot extremize F. 



3.3 Upgrade of the modulated RSB-schemes 

While the modulated i?S'i?-schemes are initially designed 
for applications to glassy ferrimagnetic systems, one should 
also consider them under a more general point of view: 
they can be used in systems without sub-lattices and even 
without reference to antiferro- or ferrimagnetic order. One 
can upgrade the schemes as an alternative of the Parisi- 
scheme provided the traces of higher powers Q'^,k > 3, 
also yield a non-divergent replica-limit. This separate point 
as well as higher order RSB is not the issue of the present 
article, where only trQ^ is needed for the KS-mode\. 



We also define a modulated scheme MMS which incor- 
porates a fc-fold repeated 5M5-structure (of size ni x ni) 
along the diagonal of an n x n host matrix, for example 
(4,6,24,3/4) where fc = 2. The ^'MS'-matrix size is cho- 
sen as a variational parameter, kept finite while the replica 
limit n — > is taken, and finally varied to extremize F. 
Altogether mAymB, and ni are variational parameters, 
which determine the embedding factor j{mA,'m-B,n 
ni) according to Eq.®, or rriA, ttib and 7 are varied and 
their selfconsistent solutions yield ^1(7) = mArriB / {rriA — 
niB + rriAj)- The free energy is obtained by means of 
liuin^Q S{n,ni)/n, where 

S'(n,rii) = E:{'j{ni))n/ni + n{n - ni){q2A + q2Bf (5) 

or S'(n, 7) = =■(7)71/711(7) + n{n - ni{'y)){q2A + 925 )^ in 
case 7 is varied. An example for a three-fold embedding is 
given by FigI21 Since relation (O is part of the definition of 
the MMS and since the replica limit maps the Parisi-type 
parameters rriA, m^, ni from [1, 00] into the interval [0, 1], 



4 Application to a layer model with a 
confined magnetic field 

In order to examine the specific features of the modulated 
S'MS'-scheme in comparison with the unmodulated one, 
we performed a detailed analysis for the 5- model (1) 
in a spatially confined magnetic field Hp. For an equal 
number of aN of A- and B-spins we chose h 
h — for the remaining (1 — a)N spins. 
The free energy of this model realization in 1-step RSB 
can be decomposed into three parts 



Hp and 



F = + aFi{Hp) + (1 - a)FiiHp = 0), 



(6) 



where 



Fo = -JafMAMB 



AT 



lim ■ 



qin) 



(7) 
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and 

F^{H,) = - J2 — r In r cosh™»{^i/.(i/p)} 

(8) 

with = /^^ dze'^"/yV2^, and Q = {Qa, Qb, Qs)- 

The effective field iJ^ on sublattice k depends on the order 
parameters of the complementary sublattice R. It is given 
in terms of 1. the confined polaron field i/^, 2. in terms of 
the magnetization field of the complementary sublattice k, 
and 3. spin fields representing the spin glass field Zk,K (in 
the 1 i?5' -B-discretized approximation), by the expression 

(9) 

The motivation for choosing this specification of the KS- 
model is essentially twofold: 

1) by scanning the full range < a < 1 we found 
that the phase diagram is not only marked by a continu- 
ous spin glass - ferrimagnet transition. At slightly higher 
ratios Jaf/J, a subsequent small flop transition from fer- 
rimagnet to an (what might be called) antiferrimagnetic 
phase occurs, which provides an ideal test-ground for the 
SMS'-scheme: 

the A <-> i3-symmetries are strongly broken, since solu- 
tions are far away from either Ma = Mb or Ma — —Mb- 
The main features of the phase diagram for all a at se- 
lected characteristic polaron fields Hp and temperatures 
is analyzed below in|21 

2) CdTe/Cdi-xMrixTe-layers are well described by 
the present model, where Hp represents a confined polaron 
field created by polarized exciton-hole spins being local- 
ized at the interface. The penetration depth of the hole- 
wavefunction defines the portion a of the CdMnTe-layer 
which is exposed to the field Hp. The magnetic CdMnTe- 
layer employs Villain- Ising pseudo-spins S ^Jj represent- 
ing tetrahedra of manganese Heisenberg-spins which re- 
tain only two orientational degrees of freedom. All pseu- 
dospins of the magnetic layer are then coupled by a long- 
range partially frustrated interaction with antiferromag- 
netic mean value. This model (with a = 0.5) provided 
optimal fits for experiments in the spin glass regime at 
X = 33% 112] . Increasing Mn-concentration x enhances the 
antiferromagnetic bias and eventually leads to a transition 
from spin glass to antiferromagnetic or ferrimagnetic order 
(in a homogeneous polaron- or external field) at a criti- 
cal concentration Xc- Scaling and numerical analysis on 
the basis of anisotropic Heisenberg models were also pro- 
vided [TMTl] . The virtue of the Pseudo-Ising concept lies 
in the smaller lower critical dimension when compared to 
Heisenberg systems. 

We explore at T = the difference between the SMS- 
scheme and the unmodulated one. For its demonstration 
we choose a polaron-field strength Hp — 4 J and a confine- 
ment-fraction a = 0.5 (half-penetrated layer). Our 1-step 
i?5'i?-results are obtained by solving up to seven coupled 
selfconsistent integral equations (S'M^'-scheme) which ex- 
tremize F. Thanks to the T 0-limit one integration can 
be solved exactly, simplifying the selfconsistent set of ten 




Fig. 3. Solutions for T = and a = 0.5 in 1-step unmodu- 
lated RSB: Parisi parameter aA = as, clk = mK./T, sublattice 
magnetizations Ma , Mb , spin glass order parameters q2A , <12B , 
and 5q2 = q2B — <12A, as a function of r = Jaf/J and finite 
confined field Hp = 4J. Dash-dotted lines locate the Ist-order 
transition, dashed lines show unstable solutions. 




Fig. 4. Corresponding results in the modulated J?S'_B-scheme 
SMS showing {aA,aB), {Ma, Mb), {q2A,q2B), and the SMS- 
embedding parameter 7. It is seen that and q^ show only 
very weak sublattice splitting. 



coupled double-integral equations for finite T considerably 
(while qiA = qiB = 1 at T = 0, both gi„(r), k = A,B, 
and the hole polarization too must be determined selfcon- 
sistently for finite T). The T = 0-resuhs of Figs|3|3|show 
a continuous SG-ferrimagnetic transition with order pa- 
rameter Ma — Mb to occur at r = Jaf /J ~ 3.25, followed 
by a discontinuous transition to antiferrimagnetic order 
near r « 4.02. As Figure |31 shows, q2A and q2B undergo 
large jumps and become (almost) interchanged at the dis- 
continuous transition in the framework of the standard 
unmodulated scheme, ie under the condition a a — clb, 
Ok = lim^^gm^/T. 

In the modulated SMS-scheme, FigQ] the selfconsis- 
tent solutions for aA,aB and q2A,q2B are different but 
show only small sublattice splitting; at the transition — 
a A changes sign, in contrast to qB — qA- 

The sublattice effective fields of EqEl help to explain 
the origin of the discontinuous transition: the competition 
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Fig. 5. a) Energies F{T = 0) of modulated IRSB- (arrow), 
unmodulated IRSB- (P) , and TJS-scheme (RS) with ferromag- 
netic /JS-solution (dashed), shown as a function of r = Jaf/ J 
in the double transition region, b) energy difference between 
modulated- and unmodulated scheme; inset (top right) shows 
susceptibility JiSB-contributions Sx- 



between JafM^ and the polaron field Hp leads (in the 
clean limit) eventually to a total spin reversal on one sub- 
lattice. Random magnetic order would smear the jump in 
any homogeneous field but the discontinuity reappears due 
to the competition between antiferromagnetic {{A, B}- 
symmetric) order, preferred in the Hp — 0-region, and 
a strongly (A, i3)-asymmetric ferrimagnetic order for suf- 
ficiently large Hp/ J. 

We finally compare our li?S'i?-rcsults for the free ener- 
gies in the double transition regime of FigsEEl Fig[51pro- 
vides evidence for the discontinuous transition to involve a 
transition from unmodulated to modulated RSB. As dis- 
cussed in higher energies correspond to improved solu- 
tions (unless identical stability criteria are met). FigEt) 
shows that the energy for the SM S-schexne is higher for 
Jaf > Jaf ~ 4.02 (lower if < holds). This crossing of 
energies at the discontinuous transition and the charac- 
ter of the i?S'i?-schemes suggests that modulated RSB 
governs the Jaf > J^j-regime. The iii'-solution is low- 
est but unstable everywhere. Despite small energy separa- 
tion, equilibrium- and nonequilibrium linear susceptibility 
shows large differences which depend strongly on the type 
of RS B-scheme (inset of FigEJ. Further observable mani- 
festations and thermal behavior remain to be considered. 



discontinuity reappears due to the competition between 
antiferromagnetic {{A, i?)-symmetric) order, preferred in 
the Hp — 0-region, and a strongly {A, i3)-asymmetric fer- 
rimagnetic order for sufhciently large Hp ^ 0. For half- 
penetrated layer (a = 0.5) the continuous transition with 
order parameter Ma — Mb, caused by the competition 
between spin glass and antiferromagnetic order, exists at 
Jaf = Jaf ~ 3.25J. To understand the a = 0.5-scenario 
in the context of all < a < 1 we found sufficient to ana- 
lyze the stability limits at T = in an i?S'-approximation. 
Introducing the definitions 



^^r M^,hp = Hp/J 



(10) 



one may cast the equation of state into the compact nested 
form 

f = uiuiO + ?? - hp)) + 77 m(m(0 +t]u{S,- hp) + hp), 

(11) 

where 

u{0 EE (1 - a)r erfi^/V2),r, = a/(l - a). (12) 
The stability limits are obtained under the constraint 
dr{M^) ^ (Pr{M, 



dM, 



0, ■ 



dM2 



^0. 



(13) 



in terms of the inverted solution t^AIk) displaying the in- 
teraction ratio r = Jaf /J as a function of the sublattice 
magnetizations M^. The constrained solutions M^ are de- 
rived by scanning all a and shown in Fig|Slfor typical val- 
ues of Hp (stability limits in terms of Jaf{oL) are omitted 
for brevity) . At a fixed a, the existence of one or three solu- 
tions implies a single continuous or a single discontinuous 
SG-ferrimagnetic transition respectively, while five solu- 
tions are necessary to obtain a double transition regime 
(with a magnetization-curve M{r = Jaf /J) shaped as in 
Figs l3l4(l . Corrections in the effective field H originating 
from intra-sublattice interactions do not change qualita- 
tively the results. The right part of FiglHI shows a special 
point which emerges for large Hp near a « 0.5 and small 

The Korenblit-Shender model in a field h 
retrieved at a = 1. 



6 Conclusions and open ends 



Hp can be 



5 Confined field-fraction effect on the 
multiplicity of phase transitions 

The sublattice effective fields of EqEl help to explain the 
origin of the discontinuous transition: the competition be- 
tween JafMft and the polaron field leads (in the clean 
limit) to a total spin reversal on one sublattice. Random 
magnetic order would smear the jump in any homogeneous 
field, as can be deduced from Fig E) at a = 1, but the 



In this article hybrid i?5i3-schemes which combine tree- 
and wave-like structures were defined. A variational em- 
bedding factor of Parisi block matrices appeared as a char- 
acteristic ingredient of the new schemes. The application 
showed that the S'M5'-modulated phase is preferred be- 
yond a critical Jaf where a type of antiferrimagnetic order 
prevails. The MMiS-scheme (in contrast to SMS) allows 
a continuous crossover to unmodulated RSB and should 
be analyzed as a candidate for glassy antiferromagnets. 
A generalization of our trace-formula to all powers of the 
order parameter matrix can create an extension of the 
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Fig. 6. Magnetizations M'^^^ g obeying the constraint 
dAI^/dr — oo, r = Jaf/J, are shown in f?S'-approximation 
for T = as a function of the polaron fraction a for field 
strengths Hp = 4, 6; the vertical line indicates the position of 
the double transition analyzed in lJ?5B-schemes in Figs l3l4l 

Parisi scheme, modelling perhaps excited states in short- 
range spin glasses. Higher orders of both types of RSB 
are currently under study, using for example techniques 
applied to unmodulated RSB in Ref.|15|. 
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